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A local order parameter which is important in the analysis of phase transitions in frustrated 
combinatorial problems is the probability that a node is frozen in a particular state. There is a 
percolative transition when an infinite connected cluster of these frozen nodes emerges. In this 
contribution, we develop theories based on this percolation process and discuss its relation to con- 
ventional connectivity percolation. The emergence of frozen order may also be considered to be a 
form of constraint percolation (CP) which enables us to draw analogies with rigidity percolation 
and its associated matching problems. We show that very simple CP processes on Bethe lattices 
lead to the replica symmetric equations for KSAT, coloring and the Viana-Bray model. 



PACS numbers: 89.75.-k, 05.50. +q,75.10.Nr 



INTRODUCTION 



There is intense interest in the relations between sta- 
tistical physics and computational complexity, from both 
the computer science and physics communities This 
activity has resulted in the application of physics meth- 
ods to computer science j3| and clever extensions of com- 
puter science methods to glassy problems 1^. One fasci- 
nating result which has emerged from these studies is the 
existence of phase transitions in computational complex- 
ity. These phase transitions are continuous in some cases 
and are discontinuous in others 0. More recently the 
A:— connectivity and A:— core problems have attracted 
interest, for example in designing redundant networks 0- 
A— connectivity is the generalisation of the conventional 
connectivity percolation problem to the requirement of 
A:— fold connectivity. That is, a graph is fc-connected if 
for each pair of vertices in the graph there exist at least k 
mutally independent paths connecting them. The A— core 
of a graph is the largest subgraph with minimum ver- 
tex degree k. The Bethe lattice equations for the fc-core 
were actually first derived in the context of A-bootstrap 
percolation[8]. A;-bootstrap percolation is the percola- 
tion process found by recursively deleting ah nodes which 
have connectivity less than k. More recently the Bethe 
lattice fc— core equations have been used to develop the- 
ories for rigidity percolation (fll ll(ilTll |. In this paper, we 
give a brief introduction to the connectivity percolation 
and g— rigidity equations on Bethe lattices and then de- 
scribe similar percolation processes which are important 
in the Viana-Bray spin-glass model, the coloring problem 
and the K-SAT problem. For these problems we develop 
equations for the probability that a infinite frozen clus- 
ter emerges. We then show that in the simplest approxi- 
mation, this formalism reproduces the replica symmetric 
equations in a surprisingly straightforward manner. 

Frozen order is a unifying concept in the analy- 
sis of glasses and geometrically frustrated systems in 
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physics [13 and also in NP-complete problems in com- 
puter science, such as coloring^^ and K-SAT|,,^. Frozen 
long-range order is most easily understood at zero tem- 
perature. At zero temperature the paradigm geometry is 
to fix the variables on a surface of the system and then 
to test whether these frozen degrees of freedom cause 
the propagation of frozen order into the bulk of a sam- 
ple. A spin is frozen only if the spin is fixed or con- 
strained by the spin configurations of its neighbors, as we 
shall demonstrate explicitly below using the Viana-Bray 
model. Frozen order may occur even though the variables 
(e.g. the spins in a spin glass) at each vertex of a graph 
look random. Furthermore, not all of the variables in the 
system need to freeze. However for the system to be in 
the frozen ordered ground state, the frozen component 
must percolate. 

The vertex q-coloring problem is equivalent to fin ding 
the ground state of the q-state Potts antiferromagnet L14( . 
Each node of a complex graph may have any one of q col- 
ors. The objective is to find the color configuration which 
minimizes the number of edges which have the same color 
at each end. The propagation of frozen color has many 
conceptual similarities with the propagation of rigidity in 
central force networks However there is a key dif- 

ference which makes the coloring problem NP-complete 
whereas the rigidity problem is polynomial. The key dif- 
ference is that the constraints in coloring are distinguish- 
able while the constraints in rigidity percolation are not. 

Spin glasses and many frustrated antiferromagnets 
map exactly to problems in the NP-complete class[J|. 
NP-complete problems are of central interest in computer 
science (CSE)(l7l| and have motivated many attempts to 
design quantum algorithms for their efficient solution. 
The phase transitions which physicists study often corre- 
spond to a change in the computational complexity of the 
corresponding CSE problem. Since these problems are of 
enormous interest in physics, CSE and also in practical 
applications it is not surprising that there is a burgeon- 
ing of efforts to better understand the phase transition 
which occurs in NP-complete problems. 

It is necessary to consider the effects of randomness 
on physics problems as randomness is present in most 
magnetic and electronic materials. The CSE interest in 
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random instances is from a different perspective. The 
motivation is to find "typical" problem instances which 
are then used to test the algorithmic complexity of new 
algorithms. A result of broad importance is the obser- 
vation of a phase transition in computational complexity 
in random satisfiability problems Q. The key quantity is 
the ratio of the number of constraints, M, to the number 
of variables, N, and this ratio is a = M/N. For a < a* it 
is believed that random SAT problems are almost surely 
in P, while for a > a* random SAT problems are almost 
surely in NP. In addition there is a phase transition as 
measured by the number of violated clauses in the opti- 
mal solution. For a < ac, the number of violated clauses 
is of order one, while for a > ac the number of violated 
clauses is of order N. It is beheved that ac > a*. 

The physics community has applied the replica method 
to NP combinatori al p roblems with remarkable success f3, 
[llll^, 20,'2T,[2ll2l- In addition new algorithms have 
been developed based on a combination of replica sym- 
metry breaking ideas from physics and belief propagation 
ideas from the artificial intelligence community [J, [2J, [2^. 
Though the replica method is an excellent tool, it is 
quite difficult both technically and intuitively. We show 
that a simple combinatorial procedure based on perco- 
lation ideas can reproduce many of the successes of the 
replica method. The percolation process occuring at the 
phase transition can be thought of as either percolation 
of constraint or percolation of frozen order. In this con- 
tribution, we derive the repHca symmetric theories for 
K-SAT, the Viana-Bray model and coloring using perco- 
lation concepts. 

The next section of the paper gives a brief review of 
the analysis of connectivity percolation on Bethe lattices 
and random graphs, and also describes its extension to k- 
connectivity percolation. Section III describes the anal- 
ysis of the glass transition, at T = 0, in the Viana-Bray 
model. Section IV focuses on the coloring problem, while 
Section V presents an analysis of K-SAT. Section VI con- 
tains a brief summary. 



II. CONNECTIVITY AND RIGIDITY 
PERCOLATION 

Percolation on diluted Bethe lattices was analysed by 
Fisher and Essamj^, who defined the probability that 
a node is part of the infinite cluster, T. They found 
that the probability that a node is not on the infinite 
cluster, Q = 1 — T, only requires that all of its connected 
neighbors also not be part of the infinite cluster, so that, 

Q = (l-p(l-Q))" (1) 

where p is the probability that an edge is present in the 
Bethe lattice, and a = z — 1, where z is the co-ordination 
number of the Bethe lattice. Note that this expression 



may be written as, 

which is more convenient for the generalisation to rigidity 
percolation. From Eq. (1), it is easy to show that there is 
a phase transition at pc = l/a and that T (p—pc) near 
the critical threshold. The phase transition is thus con- 
tinuous with order parameter exponent one. Somewhat 
earlier, this transition was also studied in the graph the- 
ory community by Erdos and Renyi0|. They concen- 
trated on random graphs, which consist of highly diluted 
complete graphs. A complete graph is a graph where 
every node is connected to every other node. In fact 
they defined p — c/N, where c is finite and showed that 
a giant (extensive) connected cluster emerges at c = 1. 
They derived an equation for the probability that a node 
is on the giant cluster, 7. Their equation is found from 
Eq. (2), by taking the limit p = c/N, N — z 00, to 
find 7=1 — e"'^'^. Near the critical point 7 ^ 2(c— l)/c^ 
so, as expected based on the universality hypothesis, 7 
also has an order parameter exponent of one. 

Rigidity percolation on Bethe lattices, is described by 
a simple generalisation of Eq. (2). In this generalisa- 
tion, each node has g degrees of freedom. For example 
if we wish to model rigidity percolation on central force 
networks, then g = d, where d is the lattice dimension. 
In order to make a giant g-rigid cluster, we need to con- 
strain the g degrees of freedom at each node with at least 
g bonds, so we generaHse Eq. (2) to, 

Ta-J2(f)(pTaya-pT,r-' (3) 

i=g ^ ^ 

which is the simple generalisation of Eq. (2) to the re- 
quirement of at least g— neighbor connections. 

Eq. (3) was first invented in the context of a Bethe 
lattice theory for Bootstrap percolationQ and has been 
used more recently to develop a Bethe lattice theory for 
rigidity percolation 0, 0|. In the random graph 
limit, Eq. (3) reduces to, 

7, = l-e--g(^ (4) 

When 5 = 1 this gives the Erdos-Renyi result for the 
emergence of a giant cluster in random graphs, while for 
g > 1 , there is a discontinuous onset of a finite solution 
at a sharp threshold Cg\^. Numerical solution of Eq. (4) 
indicates that for g = 2, C2 = 3.3510(1). This value has 
also been found in a recent mathematical analysis ;6] of 
the threshold for the emergence of the giant 3-core on 
random graphs. The k-core problem is equivalent to the 
k-bootstrap percolation problem. However the k+l-core 
is, in general different than the k-rigidity problem, and 
even on Bethe lattices and random graphs there are some 
important differences. 
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The most important diflerence is that for q- rigidity, the 
finite solution Tg is metastable for a range of c > Cg 0, ■ 
The true rigidity transition actually sets in at Cr > c„ and 
is identified using constraint counting arguments flU,'^. 
Nevertheless the probability of being on the infinite rigid 
cluster is correctly found from Eg (4) pr ovided c > c^, 
where Cr is the rigidity threshold [lO, As we shall 

see below the analogous theories for glassy combinatorial 
problems, in particular the Viana Bray model, K-SAT 
and q-coloring, provide solutions at the level of the replica 
symmetric theory. Moreover, as will be described else- 
where, the methodology we introduce here can be used 
to develop simple and accurate recursive algorithms for 
these glassy problems on general graphs. In the case of 
first order transitions, as occurs in q-coloring (with g > 3) 
and for K-SAT {K > 3), the transition point we find be- 
low may mark the onset of metastability. In order to 
find the true threshold we need to find the ground state 
energy from the order parameter, in a manner similar to 
that used in rigidity percolation. This is non trivial and 
will be elucidated elsewhere. 



III. VIANA-BRAY MODEL 

We first analyse the onset of frozen order in the Viana- 
Bray{VB) spin-glass model jS^I, which provides a basic 
model for disordered and frustrated magnets, such as 
EuxSri-xS^^. The Hamiltonian for the VB model is. 



H — ^ JijSiSj 



(5) 



where Si = ±1. The exchange constants Jij are ran- 
domly drawn from the distribution, 

DpiJv) ^p[^S{J,j+J) + ls{J,,-J)] + {l-p)S{J.,,), (6) 

As above we focus on the random graph limit p = c/N. 
We introduce the following probabilities: 

P = probability a site is frozen in the up state 
M = probability a site is frozen in the down state 
D = probability a site is degenerate 
In the absence of an applied field and within a symmetric 
assumption, P = M and D = 1 — 2M. We then need con- 
sider only one of these probabilities. However for clarity 
and for ease of generaHsation, we continue to include M 
and P separately. In terms of these order parameters, 
the magnetisation is given by, m — \P ~ M\ and the spin 
glass order parameter is, q — P + M. The recurrence 
formula for P, using p = c/N is. 
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This is understood as follows. If a bond connects a site 
at the lower level to a site at the upper level then the 
site at the upper level wants to be frozen up: if the con- 
necting bond is ferromagnetic and the lower level spin is 
frozen up; or if the connecting bond is antiferromagnetic 
and the lower level spin is frozen down. This event has 
probability, cP/2N + cM/2N. Similarly, the probabil- 
ity that a spin at the upper level of the bond wants to 
be frozen down (negative) is given by, cM/2N + cP/2N . 
The newly added spin at the upper level is frozen up if 
there is a larger number of connections from the upper 
to the lower level which prefer the frozen up state. The 
sum in Eq. (7) is thus restricted to events of this sort. 
The event (1 - c{P + M)/N) is the probability that a 
site at the lower level in the tree is either degenerate or 
disconnected from the newly added site. In the large N 
limit, Eq. (7) reduces to. 
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-^'Hcq) (8) 



where we have used the fact that we are considering a 
case where the magnetisation m = 0. In that case, 
M — P — q/2, where q is the spin glass order param- 
eter. Iq is the spherical Bessel function of zeroth order. 
The result (8) has been found before within the replica 
symmetric solution to the Viana-Bray (VB) model(see 
Eq. (15) of Q). Thus symmetric constraint percolation 
(CP) in the VB model is equivalent to the ground state 
spin glass transition as found within the replica symmet- 
ric approach. The CP approach is attractive because 
is it is simple, it avoids the mathematical difiiculties of 
the replica method and it is physically transparent. The 
construction we have used makes it clear that simple con- 
nectivity is sufficient to ensure propagation of spin glass 
order in the VB model. Constraint percolation occurs 
at c = 1 and the order parameter approaches zero as 
q 3^(c— 1), so the CP transition in this case is contin- 
uous, with the same exponent as the Erdos-Renyi tran- 
sition. 



IV. COLORING 

Now we turn to the coloring problem. Our analysis 
centers on the probability Fi {I — 1,2, ..g), which is the 
probability that a site is frozen in color I. The probability 
Fi is given by the recursion relation. 
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1=2 



(9) 
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This formula is understood as follows. In order for a site 
to be frozen in the color "1", all of the other q — 1 col- 
ors must appear, and be frozen, on one of the connected 
neighbor sites. In addition the frozen color "1" must oc- 
cur, on these neighbor sites, a strictly smaller number of 
times than any other frozen color. The probability that 
a neighbor site is connected and frozen in color "1" is 
pFi. This event occurs s times. We thus have a term 
ipFiY for the color "1". A similar term applies for each 
of the other required q — 1 frozen neighbor colors, with 
each of them occuring ki times. We must also allow for 
the possibility of events which are not of the type pFi, 
which leads to the term (1 — pJ2 ^i)''''^^ ■ This proba- 
bility is summed from to infinity as it does not have 
to exist in a configuration in order to ensure that Fi be 
finite. Note however that (1 —pj^^i) by far the most 
likely event in the random graph limit, where p c/N. 
All of these probabilities are exclusive and independent. 
We must also allow for all ways of arranging this set of 
q+1 exclusive events amongst the a possible connections 
between our newly added site and the sites at the lower 
level in the tree. This leads to the multinomial factor. 
An equation like (9) occurs for each of the q colors which 
are allowed. If we assume that all colors have the same 
probability (which is natural provided there are no sym- 
metry breaking terms), then Fi = F2 = Fi = F/q. Using 
this, and taking the random graph limit yields, 

s=0 ^ fc=s+l ^ 

This equation is valid for arbitrary q provided q/N 0. 

For g = 2, we assume that F is continuous near the 
percolation threshold and expand this expression in pow- 
ers of F which yields, 

F^cF-\{cFf + 0{{cFf) (11) 

This has the solution, 

c>l (12) 

This is, other than a prefactor of 4/3c^ instead of 2/c^, 
the same as the critical behavior of the giant cluster prob- 
ability in random graphs [13, E^l- For c well away from 
the transition, we solve Eq. (10) numerically. The s and 
k sums are rapidly convergent and for the c range near 
criticality, only a few terms are required for high accu- 
racy results. From the solution for F we obtain all of 
the results of interest and they are presented in Fig. 1. 
The continuous behavior of 2-coloring near threshold is 
evident from these data. 

For q — i, sjn attempt to find a continuous transi- 
tion by expanding in powers of F fails. Numerical solu- 
tion of Eq. (10) is presented in Fig. 2 where it is seen 
that there is a jump discontinuity in the infinite frozen 
cluster probability at a sharp threshold. We find that 
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Figure 1: The coloring order parameters for q — 2. The 
lower two curves are the probability that a site is frozen and 
colorable, G (the s = term in Eq. (10)), and the probability 
that a site is frozen and frustrated, H (the s > 1 term in 
Eq. (10)). The top curve is the probability that a site has a 
frozen color F = G + H , which is found by solving Eq. (10) 
with q = 2. 
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Figure 2: The coloring order parameters for q = 3. The 
lower two curves are the probability that a site is frozen and 
colorable, G (the s = term in Eq. (10)), and the probability 
that a site is frozen and frustrated, H (the s > 1 term in 
Eq. (10)). The top curve is the probability that a site has a 
frozen color F = G + H , which is found by solving Eq. (10) 
with q = 3. 



c* = 5.14(1) and that the jump in the order parame- 
ter is AFc = 0.365(1). We thus find that the coloring 
transition for g = 3 is first order as has been found in nu- 
merical simulations [T^ on random graphs. Our coloring 
threshold is consistent with a recent replica symmetric 
numerical calculation, which yielded c* « 5.l[23) but 
is significantly higher than that found in the simulation 
work of Culberson and Gent[3| where c* w 4.5 — 4.7 or 
in the numerical work on survey propagationfj^ . which 
yields c, « 4.42. Nevertheless the nature of the transi- 
tion is correctly captured by the simple CP theory. It is 
also important to note that the solution found here may 
also be metastable for a range of c, as was found in the 
rigidity case^^- The onset of metastability is an impor- 
tant threshold from the point of algorithmic efficiency, 
as it marks the onset of glassy relaxation dynamics. The 
coloring theory developed above can be formulated in a 
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very similar way to the formulation of the propagation 
of the k-core. However there is a critical difference. The 
constraints in the coloring theory have to be treated as 
distinguishable, while the constraints in the k-core calcu- 
lation are indistinguishable. 

V. K-SAT 

The satisfiability problems we consider ask the ques- 
tion: Given a set of binary variables, % = 0, 1 or equiva- 
lently Zi = True orFalse, is it possible to satisfy a spec- 
ified set of constraints on these variables? In the K-SAT 
CcLSG, cL typical constraint is of the form, 

(ziAZjAzk) (13) 

where A is the logical OR operation and the overline in- 
dicates a negated variable. This logical clause is satisfied 
(SAT) if any one of the variables in the clause is SAT. 
The variables Zi and Zk are SAT if they are true (T), 
which we take to he Zi = z^ = I, while the variable Ij is 
SAT when Zj is false (F), which corresponds to Zj = 0. 
We shall also fix the number of variables in each clause 
to be K, which is the K-SAT problem. In these SAT 
problems we shall randomly choose a set of M clauses 
like that in Eq. (13)and try to find the assignment of 
the binary variables which minimizes the number of vi- 
olated clauses. Each variable appearing in a clause is 
negated with probability 1/2 and the number of vari- 
ables is N. The key ratio is a = M/N. We would like 
to find the threshold for constraint percolation. That is, 
what is the threshold for the appearance of a giant clus- 
ter of clauses where each clause is completely specified 
or "frozen". These completely specified clauses cannot be 
altered without increasing the total number of violated 
clauses, so that they are non-degenerate. There are three 
types of clauses in an optimal assignment of a formula: 
(i) Clauses that are SAT but are degenerate; (n) Clauses 
that are SAT but are frozen; (in) Clauses that are UN- 
SAT but arc degenerate. Only type (ii) clauses propagate 
constraint. We seek a formula for the probability, V, that 
a variable is frozen and the probability, F, that a clause 
is frozen and SAT. 

We make a tree construction of the factor graph for 
the K-SAT problem (see Fig. 3). The probability that 
a variable is frozen and part of the giant frozen cluster 
is V and the probability that a clause is frozen and part 
of the giant frozen cluster is F. The branching of the 
variable nodes has maximum co-ordination M, but the 
probability that a link actually exists between a node and 
clause is p = K/N. We start by assuming that a variable 
is frozen at level 1 (sec Fig. 3) and then determine the 
consequences of this assumption at levels 2 and 3. 

The probability that a clause is frozen, F, at level 2, 
given the probability, V, that a variable is frozen at level 
1 is given by, 

F = {^f-K (14) 




V = probaUlit)' a mriabk node is frozen 
F - probability a clause node is from 

Figure 3: The factor graph used to construct the recurrence 
relations. The circles denote variable nodes, while the square 
nodes are the clause nodes. V is the probability that a vari- 
able node is frozen, while F is the probability that a clause 
node is frozen (see the text). We assume that a variable at 
level 1 is frozen and find the probability that a variable at level 
3 is frozen. The clause nodes have co-ordination K, while the 
variable nodes have co-ordination M 



This equation is understood as follows. In order for a 
clause at level 2 to be frozen by the variables at level 
1, all of the level 1 variables to which it is connected 
must be frozen and in conflict with the assignment in the 
clause. This imposes a fixed assignment on the variable 
3. This is the only configuration of variables at level 1 
which propagates constraint through a clause to level 3. 
Now we must consider the cummulative effect of all of 
the clauses which are connected to the variable at level 
3. There arc M — 1 such clauses of which a fraction F 
propagate constraint (are frozen) according to the mech- 
anism of the previous paragraph. Some of these frozen 
clauses propagate the requirement x and others propa- 
gate the requirement x. The variable at level 3 then has 
three possible states, P = positive, = negative and 
£)= degenerate. The state of the level 3 variable is de- 
generate if the number of constrained connections which 
favor the positive state (.x) is the same as the number 
of connections which favor the negative state (x). The 
probability this variable is frozen (ie. either negated or 
not) is V = P + N = 1 — D as we are considering the case 
where the probability that a variable is negated is 1/2. 
It is straightforward to generafise to the case of unequal 
probabilities. The probability that the node at level 3 is 
degenerate is then, 

D = Y — (^)2'=(1 - pF)^-2'= (15) 

^ (fc!)2(M-2fc)!^ 2 ^ 

Where we have used the fact that the probability that a 
connection occurs between a variable node and a clause 
node is p = K/N. Eq. (15) is understood as follows. 
The probability that a clause at level 2 is frozen and con- 
nected (ie it propagates constraint), and it requires the 
variable at level 3 to be a; is pF/2. The probability that 
this clause propagates constraint and requires the vari- 
able at level 3 to be x is also pF/2. The variable at 
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Figure 4: The probability that a clause is frozen, F, as a Figure 5: The probability that a clause is frozen, F, as a 
function of a = a, for 2-SAT. function of a = a, for the 3-SAT problem 



level 3 is degenerate if these two events occur an equal 
number of times, hence the term . The combi- 

natorial factor gives all ways of arranging these events, 
taking into account that the x and x events are distinct. 
In the thermodynamic limit, using pM — aK , we find. 



D = 



-aKF \ " 



1 ,aKF 



2fe 



(16) 



This provides the reccurence formula for V =\ — D which 
may be written in the form. 



y = 1 - e 



-aKF 



h{aKF) 



(17) 



where /q is the spherical Bessel function of zeroth order. 
Note that /o(0) = 1- For completeness, we note that the 
probability that the new variable is frozen in the positive 
(not negated) state is, 



M M 



^ ^ ^ (fc!)2(Af-fc-0!' 2 



k=0 l=k+l 



(18) 

The probability that the variable is frozen in the A'^ state 
is the same as P for the case we are considering, where 
the variables have equal probability of being negated and 
not negated. 

Equations (14) and (18) provide the self consistent 
theory for the onset of a giant constrained cluster in 
K-SAT. We now analyse this theory for the two typical 
cases. 

The 2-Sat case (K^2) 

In this case Eq. (14) is F = V/2. Expanding Eq. (17) 
in powers of F, we then have, 

F = i[l - (1 - 2aF + + ..)(1 + a^F^))] (19) 

This has the trivial solution a = 1. It also has the non- 
trivial solution 



(20) 



Thus the random 2-Sat giant cluster emerges smoothly 
at a = 1. Numerical calculation of F from Eqs. (14) 



and (17) is presented in Fig. 4. 

The K > 3-Sat case 

In these cases, Eq. (14) with (17) do not have a non- 
trivial solution with a smooth behavior near a criti- 
cal point. However they do have a non-trivial solution 
which has a discontinuous onset at a threshold value, 
ac{K). This solution is found by iteration of Eq. (14) 
and Eq. (17) and the results are presented in Fig. 5. We 
find that although the emergence of the giant cluster is 
discontinous, for any K > 2 the size of the first order 
jump decreases quite rapidly with increasing K. This 
indicates that the K-SAT transition is weakly first order 
and that an analytic analysis at large K is possible. The 
3-SAT critical value which we find, etc (3) ~ 4.6673(3), 
is consistent with the replica symmetric solution|J2j for 
the metastability point, and significantly higher than 
the numerical values for the K-SAT transition which 
lie around 4.30. The numerical results we have found 
(using Eq. (18)) for metastability point and the jump 
in F at that point are: ac(3) = 4.6673(3), 6Fc — 
0.0680(1); ac(4) = 11.833(1), 6Fc = 0.0341(3); ac{5) = 
29.91(1), SFc = 0.016(1) ; ac(6) = 64.1(1), SFc = 
0.0071(1). 



VI. SUMMARY 

We have shown that the probability that a site is on 
the infinite frozen cluster may be calculated using sim- 
ple combinatorial methods. This provides a general an- 
alytic approach to many hard combinatorial problems, 
and provides a useful complement to the replica method. 
Although we concentrated on the symmetric theory here, 
cavity methods[l| hold promise for generalising this ap- 
proach to the unsymmetric case, as will be presented else- 
where. 

The coloring transition is continuous for g = 2 and 
discontinous for 9 > 3, similarly K-SAT is continuous 
for K = 2 and discontinuous for K > 3. In contrast 
the VB model of glasses has a continuous phase tran- 
sition. As found in the rigidity percolation problem^, 
processes which require more than 2-connectivity in or- 
der to propagate constraint have a tendency toward first 
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order transitions. However a counter example is rigid- 
ity percolation on triangular lattices, where the rigidity 
transition is continuous|2d|- It thus seems a difficult task 
to determine the conditions which produce continuous as 
opposed to discontinuous percolation transitions in com- 
plex combinatorial problems. 



discussions with Radu Cojocaru, Remi Monasson and 
Bart Selman. PMD thanks Chris Farrow for a careful 
reading of the manuscript. 
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